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ABSTRACT 


Most analytic and computer simulation models of the Navy's 
CONREP phase of underway replenishment operations have assumed expo- 
nentially distributed service times for two reasons: (1) lack of 
better estimate of the true distribution, and (2) such an assumption 
leads to more tractable mathematical solutions. 

By analyzing several different combinations of replenishment 
vessels and combatant ships, it is demonstrated through goodness-of- 
fit tests that gamma distributions whose parameters can be estimated 
from actual operational data are more precise estimates of the actual 
underlying CONREP service time distributions. Furthermore, it is shown 
that the distribution can be made Erlang by minor adjustments to the 
parameter estimates of the gamma distribution. Such a procedure might 
be desirable for analytical models employing Laplace transforms. 

The data is subjected to linear regression techniques in an 
effort to develop meaningful and accurate functional relationships 
between service time and customer needs. The results indicate that 
the standard error of the estimating relationships would be too large 
to be of any practical use in planning an underway replenishment 


operation. 








TABLE OF CONTENTS 


Section 


I INTRODUCTION 


fi, “OBJEGIIVES 


Ili. DISCUSSION OF DATA SOURCES 


IV. DESCRIPTION OF DATA ANALYSIS TECHNIQUES UTILIZED 


V. RESULTS OF DATA ANALYSIS 


A. Determination of CONREP Service Time Distributions 
and Their Parameters 


1. AO Serving CVA/CLG/DD Type Customers 
2. AE Serving CVA/CLG/DD Type Customers 
3. AOE Serving CVA/CLG/DD Type Customers 
4. FS Serving CVA/CLG/DD Type Customers 
5. AF Serving CVA/CLG/DD Type Customers 
6. AKS Serving CVA/CLG/DD Type Customers 
B. Application of Linear Regression Analysis to 
CONREP Data 
VI. CONCLUSIONS 


BIBLIOGRAPHY 


APPENDIX 
A The Gamma Distribution 
B Computer Subroutine GAMDIS 
6 Scatter Diagrams of CONREP Data 


28 


33 


34 


35 
oo 
4] 





Table 


Il 


IV 


V1 


Vii 


VIII 


[x 


LIST OF TABLES 


Statistics on Determination of Estimates of Parameters 
of Gamma Distribution for Pilot Cases 


Statistics on Determination of Estimates of Parameters 
of Gamma Distribution for Pilot Cases, r Converted to 
Nearest Positive Integer 


Statistics on AO Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Statistics on AE Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Statistics on AOE Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Statistics on AFS Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Statistics on AF Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Statistics on AKS Serving Situation, CVA/CVS, 
CLG/DLG/DDG, and DD Type Customers 


Results of Linear Regression Analysis on AO and AE 
Serving Situations 


Page 


19 


20 


25 


Zo 


26 


26 


2/ 


2/ 


30 





Figure 


WE 


i. 


1S. 


14. 


io). 


HO. 


LIST OF ILLUSTRATIONS 


Frequency Histogram of CONREP Data for AE-DD Situation 


Gamma Probability Distribution in the Special Case 
Where A = | 


BIMED 02D Scatter Diagram for AO-CVA/CVS Situation 
(Day), NFSO vs. Service Tile 


BIMED 02D Scatter Diagram for AQ-CVA/CVS Situation 
(Day), JP-5 vs. Service Time 


BIMED 02D Scatter Diagram for AQ-CVA/CVS Situation 
(Day), AVGAS vs. Service Time 


BIMED 02D Scatter Diagram for AQ-CVA/CVS Situation 
(Night), NSFO vs. Service Time 


BIMED 02D Scatter Diagram for AO-CVA/CVS Situation 
(Night), JP-5 vs. Service Time 


BIMED 02D Scatter Diagram for AO-CVA/CVS Situation 
(Night), AVGAS vs. Service Time 


BIMED 02D Scatter Diagram for AQ-CLG/DLG Situation 
(Day), NSFO vs. Service Time 


BIMED 02D Scatter Diagram for AO-CLG/DLG Situation 
(Night), NSFO vs. Service Time 


BIMED 02D Scatter Diagram for AO-DD Situation (Day), 
NSFO vs. Service Time 


BIMED: 02D Scatter Diagram for AO-DD Situation (Night), 
NSFO vs. Service Time 


BIMED 02D Scatter Diagram for AE-CVA/CVS Situation 
(Day), AE Material vs. Service Time 


BIMED 02D Scatter Diagram for AE-CVA/CVS Situation 
(Night), AE Material vs. Service Time 


BIMED 02D Scatter Diagram for AE-CLG/DLG Situation 
(Day), AE Material vs. Service Time 


BIMED 02D Scatter Diagram for AE-CLG/DLG Situation 
(Night), AE Material vs. Service Time 


Page 
17 


38 


42 


43 


44 


45 


46 


47 


48 


49 


50 


2 | 


52 


518) 


54 


Sho 


Figure Page 


ize BIMED 02D Scatter Diagram for AE-DD Situation (Day), 
AE Material vs. Service Time 56 


18. BIMED 02D Scatter Diagram for AE-DD Situation (Night), 
AE Material vs. Service Time 57 


ACKNOWLEDGEMENT 


The writer is deeply indebted to Professor Paul R. Milch 
for his guidance and encouragement through the course of this 


investigation. 





I. INTRODUCTION 


The United States Navy recognizes that one of the most prac- 
tical ways to increase the effectiveness of our naval forces in main- 
taining Superiority at sea is through the application of Underway 
Replenishment (UNREP). The purpose of these operations’ is to improve 
unit endurance by supplying needed logistic items to fleet units with 
minimum increase in vulnerability and with minimum diversion from 
performance of combat missions. The resulting ability of ships to 
remain on or near station is equivalent to adding additional combatants. 

The importance of UNREP has been reemphasized within the 
Department of the Navy in that UNREP responsibility is now a 
NAVSHIPSYSCOM designated Project Office (PMS-90) [1]. NAVSHIPS 
Instruction 5430.73 of 25 August 1967 is the implementing document. 
Prior to this the responsibility for UNREP was solely a technical one 
and located in various BuShips branches. — 

There are basically three ways in which an UNREP may be accom- 
plished. The primary method used is called Connected Replenishment 
(CONREP). This operation consists of intership horizontal transfers 
via rigs connecting the supply ships and the supported units. An 
increasingly important method, which involves the use of UH-34 or 
UH-46 helicopters to deliver certain items, is that of Vertical 
Replenishment (VERTREP). It has been envisioned that VERTREP will 
Someday enable al] but non-nuclear ships to avoid the inherent vul- 
nerability imposed during CONREP. The third method entails a com- 
bination of CONREP and VERTREP. 

Considerable effort has and is being expended to reduce the 


total time to conduct an UNREP. These efforts include extensive 





training, improved delivery techniques and equipment, and newer designs 
in supply vessels. Several attempts have been made to model the under- 
way’ replenishment operations, both from the analytic and computer 
Simulation approach. 

McCullough [2] made an analytic approximation of the replenish- 
ment process by using a multi-stage cyclic-queueing model. The model 
considered M supply ships, called stages, placed in series. These 
ships serviced N combatant units, each of which passed by the supply 
ships in succession. In order to achieve’ a tractable solution, the 
operation was permitted to repeat itself continuously in order to 
study the steady state behavior of the system. “A strong assumption 
used in this study was that of exponential services times. The most 
basic objection to this approach is that it is not really representa- 
tive of what actually happens during CONREP. First, the cyclic nature 
of this model’ is a great simplification of the more complex sequences 
possible in these operations. Secondly, the operation is finite, and 
therefore’ no steady state solution ever exists. The assumption of 
exponential service times in order to make the model’ more amenable to 
analysis is not valid, a fact that this paper will attempt to prove. 
Included’ in McCullough's work: is an incremental’ time step computer 
Simulation: model to test various’ sequences and to:make comparisons 
with the results obtained in the’analytic mode?. Here again, a basic 
assumption was’ the use of specified exponential service times for the 
Supply ships, which did not vary between different types of customers. 

Gordon’ and: Copes [3] developed a deterministic model to 
investigate the optimal scheduling of underway replenishment operations 


by considering the process as a job scheduling problem. The service 
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times of each combination of supply ship and customer were assumed 
to be deterministic and known, and general expressions were obtained 
for the total time to complete the CONREP and the total waiting time 
of the ships involved. The most serious deficiency of this approach 
is the assumption of fixed service times. It is more reasonable to 
assume that such service times are in fact random variables, and that 
questions relating to the determination of the average time needed to 
complete an underway replenishment operations must take this into 
account. 

Waggoner [4] has developed an analytic model for the case of 
two supply ships and L combatants using queueing theory concepts and 
a random walkin the plane. The distribution of the total time to 
complete the finite operation is obtained in terms of its Laplace 
transform, which although frequently uninvertable, can provide informa- 
tion regarding the mean and variance of CONREP completion times. Here 
again, nowever, exponential service times were assumed. Waggoner 
concedes the point that although the exponential distribution fre- 
quently fits many realistic queueing models, some doubt exists as to 


its applicability in the case of CONREP operations. 








ioe OBJEGIIVES 


The objectives of this paper are twofold: 

1. To determine from operational data more accurate estimates 
of the distribution of service times between various participants in 
the CONREP situation. If this can be accomplished, then more accurate 
results from future analytical and computer simulation models on CONREP 
operations can be expected. 

2. Through the use of scatter diagrams and statistical regres- 
Sion techniques, determine if functional relationships extst between 
the amount’ of a specified type of material received by the combatant 


and the service time necessary to provide this material. 


14 





TV. GENERAL’ APPROACH. TO CONREP DATA ANALYSIS 
A. DETERMINATION OF CONREP SERVICE TIME DISTRIBUTION AND PARAMETER 
ESTIMATES 

From previous experience, it was felt that. service times in a 
CONREP operation would vary in at least two ways. First, the total 
requirement and replenishment configuration would differ among com- 
batants serviced by the same replenishment ship. Secondly, variation 
would occur among servers replenishing the same type of customers, 
since delivery rates vary from server to server. 

One of the first steps taken to determine the distribution of 
service times was’ to’ look carefully at the data from several specific 
Situations, and from the results find a pattern which might make it 
possible to make“a’ statement regarding the underlying distribution of 
CONREP service times in’ a general sense. 

The three pilot cases chosen were those consisting of replenish- 
ment ship types AOE, AFS, and AE servicing DD (Destroyer) type customers. 
These cases were selected primarily due to the large number of sample 
points available, but also because it would permit an evaluation of the 
hypothesis held a priori. 

Frequency histograms were constructed from these sets of data, 
one of which is shown in figure 1. Because of the distinctive skewed 
Shapes of these graphs, it was decided to investigate the possibility 
of fitting a "gamma" distribution to each of these situations. The 
gamma distribution is discussed in detail in Appendix A. 

Since the two parameters of the gamma distribution can be 
estimated from the sample mean and variance, these statistics were 


computed (utilizing appropriate computer subroutines). With these 
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parameters estimated, another subroutine developed for this analysis 
called GAMDIS was utilized to provide the necessary computations neces- 
Sary to conduct goodness-of-fit tests with the observed data. Since 
the subroutine GAMDIS is a rather unique program involving numerical 
integration and is not to be found in the usual reference literature, 

a copy is provided in Appendix B. 

The chi-square goodness-of-fit test was used as a measure of 
the appropriateness of fit. Having the expected and observed numbers 
of units possessing service times which fall into "cells" which 
represent an exhaustive and mutually exclusive set of intervals for 
which data is held, the problem is to test the hypothesis 

H: Ps=Pa, Cine ark) 
where OF is the probability of finding a service time in the Tey cell, 
pe is the observed relative frequency of the actual data falling 
within the same intervals, and k denotes the number of cells. Ten 
minute cells were utilized, with the exception of the last cell, which 
had an interval of time ranging from the end of the final cell having 
a data point within’ it to plus infinity. The test procedure is to 
calculate the chi-square statistic 

k 


aa Z 
xe =) COp-ES)I7E, 


where 0. represents the number observed in the gs cell and E.=np. 
represents the number expected in cell i if H is true, with n equal 
to the sample size. The decision rule is 
reject H if x’ > X* (1-0) (kB) 
It is necessary to use k-3 degrees of freedom since the two population 


parameter estimates are obtained from the data. 


The results obtained employing the above procedures on the 
three pilot cases are shown in Table I. In the table, MEAN ue denotes 
average service time, VAR. is the sample variance, STD. DEV. is the 
sample standard deviation, Yr and X are estimates of the parameters of 
the gamma distribution as computed from the data, x* is the chi-square 
Statistic, and XK 95 is the tabled chi-square value with k degrees 


of freedom, and significance level of a = .05. 


Table I. 


Statistics on Determination of Estimates of Parameters for 
Gamma Distribution for Pilot Cases 


ME AN SOT 
SERVER lie VAR. DEV. Yr 


>"> 


2 2 
x Xk, 95 


COMBINED 210) 5 462.7 Za oe: al 24.5 SoZ 


AOE 48.0 308.0 lds rims mane Zine AST) 
AE 90.6 Soca oil, ae /@2 Calused, 302 


AFS on Zee evn daliae aa 8.4 Zo a0 


These results indicated that, for the particular cases chosen, suitable 
parameter estimates could be obtained from the data which enabled gamma 
distributions to be fitted to each situation. As indicated, the cal- 
culated value of the chi-square is somewhat less than the tabled value 
in each case. | 

As noted in Appendix A, the gamma distribution becomes Erlang 
if the parameter r is a positive integer. Since the Erlang is somewhat 
easier to use in analytic work involving the use of Laplace transforms, 
it was decided to adjust the value of r in each case to the nearest 
integer, and then find the value of X which minimized x2. The results 
of this procedure as applied to the three pilot cases are contained in 
Table II. 
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Table II. 


Statistics on Determination of Estimates of Parameters for 
, Gamma Distribution for Pilot Cases, 
r Converted to Nearest Positive Integer 


SERVER r N E X? 6 
COMBINED 5.0 WW 26.2 35.2 
AOE 7.0 14 22.3 26.7 
AE 4.0 .08 19.4 35.2 


AFS 12.0 nZS Gro 2580 


It is interesting to note that an improvement in X? was ob- 
tained in all but one case. For this reason, and also due to the 
advantage of working with integer values of r whenever possible, the 
practice was followed throughout the remainder of the effort associated 
with service time distribution computations. Since it was always 
necessary to compute the actual value of the Yr estimate and its 
associated X*, it was possible to note that the practice improved 
(lowered) the X? value approximately 60% of the time. It never raised 
the X* high enough to cause it to fail a goodness-of-fit solely because 
of this adjustment. 

It should be mentioned at this point that no consideration had 
yet been given to the possibility that the distribution of service time 
might also vary depending upon whether the CONREP was conducted during 
daylight or at night. Eventually surmising that this might be the case, 
the data was separated and the analysis rerun. As expected, the average 
service time for night time operations was generally higher, and adjust- 
ments were necessary in the estimated parameters r and \ in each case. 


Therefore all remaining calculations were conducted on data which had 
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been carefully separated between day and night operations. Fortunately, 

all data sources clearly indicate this important factor. 

B. LINEAR REGRESSION MODEL TO DETERMINE THE RELATIONSHIP BETWEEN CONREP 
SERVICE TIME AND AMOUNT OF MATERIAL TRANSFERRED 

Prior to undertaking this part of the data analysis, it was felt 
that nearly linear’ relationships existed between customer needs and the 
times necessary to provide them. If* these relationships could be found, 
they could be a valuable aid to those responsible for planning an UNREP. 
For example, if this technique could be utilized to determine the total 
time of an UNREP, then this: knowledge would aid planners and partici- 
pants alike since more effective use of combatants and Underway Replen- 
ishment Groups (URGs) would result. 

In order to provide a visual aid in determining what relation- 
ships actually existed between service times and requirements, scatter 
diagrams were produced using BIMED 02D [10], a statistical subroutine 
contained in the function library of the IBM 360 computer facility at 
the Naval Postgraduate School. The cases chosen as representative of 
the many possible combinations were those involving AQs and AEs as 
replenishment vessels’ and combatant types CVA, CLG, and DD. This 
selection was chosen since AOs provide an assortment of liquid fuels, 
while AEs supply only dry material, namely, ammunition and explosives. 

Before discussing the scatter diagrams further, it is neces- 
Sary to digress for the moment to describe the unique problems 
presented by the AO server situation. These supply vessels have the 
capability of supplying three different fuels, namely, Navy Standard 
Fuel Oi] (NSFO), Diesel Fuel (JP-5), and Aviation Gasoline (AVGAS). 


One customer, the conventionally powered aircraft carrier, frequently 
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takes all three types of fuel simultaneously through separate hoses. 
Other customers, say helicopter equipped, conventionally powered 
cruisers and destroyers, will take NSFO or JP-5 and perhaps AVGAS. 
As can be seen on the scatter diagrams, there is a wide range of fuel 
type and requirements’ combinations. Thus the problem of finding a 
Single, controlling independent variable for the AO server situation 
would’ appear to have no solution. The alternative is to attempt to 
find a method of predicting the response based on the magnitudes of 
a number’ of’ variables. 

The scatter diagrams for 16 situations involving the afore- 
mentioned participants are shown in figures 3 through 18 in Appendix C. 

In the special case of the AO serving CVA/CVS ships, each of 
the independent variables, barrels of NSFO, JP-5, and AVGAS, were 
plotted individually against total CONREP service time. In all other 
cases, only the dominant independent variable was plotted with time. 
All: daylight CONREP’ situations are plotted separately from night 
operations. 

A discussion of the scatter diagrams is contained in the 
next section. 

The technique utilized to discover the functional relation- 
ships between needs and service times was that of linear regression. 
In the AO-CVA/CVS cases, usually involving three independent variables, 
multiple linear regression was applied. The aim was to determine if 
a regression model of the form 

Ls =q + 8X + BX, + BX, +E 
could be calculated and applied in this case, where Lf represents the 


response dependent variable, total service time, and Xs Xo and x, 
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represent the independent variables, barrels of NSFO, JP-5, and AVGAS 
respectively. E denotes the error term in the model. 

Since the remainder of the cases involve only one independent 
variable, the attempt was to find a simple linear regression model of 
the form 

bis =O te 
where X represents barrels of NSFO for the other AO serving situations, 
and short tons of explosive material for all AE serving situations. 

BIMED O03R, a multiple regression subroutine contained in the 
function library of the IBM 360 computer facility at the Naval Post- 
graduate School, was utilized for all regression work. 


No attempt was made to fit higher order polynomials to the 


CONREP data. 


Zs 


V. RESULTS OF THE DATA ANALYSIS 
A. DETERMINATION OF CONREP SERVICE TIME DISTRIBUTIONS AND PARAMETER 
ESTIMATES 

The techniques described in Section IV for the determination 
of the parameters of the gamma distribution from the data were applied 
to a total of 29 out of a possible 36 separate cases, consisting of the 
combinations of AO, AE, AOE, AFS, AF, and AKS serving CVA/CVS, CLG/DLG/ 
DDG, and DD type customers, for day and night operations separately. 

The results of this determination are contained in Tables III 
through VIII. The following abbreviations are used in the tables. 

N = total number of sample points considered or 


available. A maximum of 30 was utilized to 
reduce size of problem. 


vm = average CONREP service time of sample. 
ons = sample standard deviation of CONREP service time. 
r = estimate of parameter r for gamma distribution 
to nearest integer. 
h = estimate of parameter X for gamma distribution. 
xe = chi-square statistic for goodness of fit. 
k = degrees of freedom. 
Pe 95> tabled value of chi-square statistic with k degrees 


of freedom and significant level a = .05. 


Statistics were not computed for situations in which eight or 
less data points were held. 

Only two cases, AO-DD (Night) and AFS-CVA/CVS (Day) failed the 
chi-square goodness-of-fit tests at the a = .05 level. Both would have 
passed at the .03 level. 

Eight out of 29 cases had estimates of the r parameter of 9 or 


greater, indicating that these cases could be estimated by a normal 
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population with parameters of mean and standard deviation indicated 


in the tables. 


prs 
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Table III. 


Statistics on AO Serving Situation 


CVA/CVS CLG/DLG/DDG 
DAY NIGHT = DAY_—=—Ss NIGHT 
30 a2 30 30 
114.3 123.4 Sonu ian 
34 .] 30.4 Zee 22m) 
11.0 16.0 720 110 
09 mules: ml fale 
20.5 15.9 202 9.9 
19 17 le 12 
S07 27.6 22.4 21.0 
Table IV. 
Statistics on AE Serving Situation 
CVA/CVS CLG/DLG/DDG 
DAY NIGHT = DAY_—sNIGHT 
30 30 30 30 
98.6 161.2 7158 1e8ig3 
e.2 vas Sie5 = zemn 
one 9.0 2249) 2.0 
703 06 FOZ 01 
ZONG 27.0 34.7 62.4 
26 27 26 5] 
3029 40.1 33.9 68.7 


DD 


DD 


NIGHT 
30 
78.3 
29.6 
7.0 
.09 
2ena™ 
17 
27.6 


NIGHT 
30 
83.3 
36.1 


Table V. 
Statistics on AQE Serving Situation 


CVA/CVS CLG/DLG/DDG DD 
DAY NIGHT DAY NIGHT DAY NIGHT 


19 5 23 3 30 11 
( i 
141.7 1 81.2 1 51.6 7089 
( | 
58.9 ra 41.1 rm 18.3 35.7 
<r <— 
SS [aan 
6.0 ia 4.0 — 8.0 4.0 
Pas = 
Lij | 
.03 O 005 oO 17 £05 
= > 
28.0 5 27.1 =) a7... 18.2 
VY VY 
= = 
31 1 17 8 13 
i] l 
45.0 27.6 5 2022 a 
Table VI. 


Statistics on AFS Serving Situation 
CVA/CVS CLG/DLG/ DDG DD 
DAY NIGHT DAY NIGHT DAY NIGHT 


13 3 14 13.30 24 
t 
145.8 1 10am3 = 7D. 52m4 58. 
61.6 = 60.0 30.3 12.2 21.3 
So 
5.0 = 3.0 5.0 18.0 7.0 
o3 © .02 07.35 m2 
3 | 
ps5" 3 73 #20 gee 4.3 
<= 
30 q 26 |Z 6 9 
( 
43.8 1 38.9 21.0 12.6 16.9 
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Table VII. 
Statistics on AF Serving Situation 


CVA/CVS CLG/DLG/DDG DD 

DAY = - NIGHT =DAY_ =—-NIGHT = DAY_—s NIGHT 
N 9 6 16 9 30 20 
T ion | 61.9 69.2 55.2 54.3 
one 45.3 : 29.2 27.7 ae 
r 6.0 = 4.0 B70. 0M 6.0 
d 05 3 06 08 18 11 
x2 25.9 i 16.7 9.1 ge aie 
k 17 7 13 1 7 9 
X95 27.6 ! 22.4 19.7. 14.1 16.9 


Table VIII. 
Statistics on AKS Serving Situation 


CVA/CVS CLG/DLG/DDG DD 
DAY NIGHT DAY NIGHT DAY NIGHT 
N 3 6 20 5 30 19 
= | 
ij 41.3 36.2 45.9 
S 
“A i ] < 
Tee a. < 13.5 = lieo — oma 
rR = a a 
r ras = 9.0 I 9.0 5.0 
“. S a a 
r “ Al oO 24 .10 
= i = 
x2 in rm 2.7 = 6.3 8.0 
= Y = 
k z ~ 6 1 6 9 
| i | 
i i i 
XV 96 : 1 12.6 1Zuom «6=Oo16.9 
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B. LINEAR REGRESSION MODEL TO DETERMINE THE RELATIONSHIP BETWEEN 
CONREP SERVICE TIME AND AMOUNT OF MATERIAL TRANSFERRED . 


With few exceptions the scatter diagrams in Appendix C tend to 
discourage the notion that meaningful, accurate linear relationships 
exist between combatant ship requirements and service time. This is 
especially true in the case of the AO servicing various units. There 
is a definite improvement in those situations involving the AE as the 
server. On several of these diagrams, it is possible to fair in a 
line which would serve as a linear model. It is questionable from 
the results’ of the scatter diagrams that the variability around a 
Single’ regression’ line could be low enough to be acceptable. In order 
to determine this’ variability°more precisety, the multiple regression 
techniques mentioned in Section IV were employed. The results of these 
efforts are shown in Table IX, for the AO and AE server situation only. 

As’ an example, this table indicates that a linear regression 
model for predicting the amount of time necessary for an AO to supply 
quantities of NSFO, JP-5, and AVGAS could be computed using the 
following formula: 

UP = 68.1 + .0050 x + .0014 x, + .0062 x. (minutes ) 
where x denotes barrels of NSFO, x, barrels of JP-5, and x, barrels 
of AVGAS required. The standard error (i.e., the standard deviation) 
of the estimate is 27.0 minutes. Three cases were selected from the 
dependent sample to illustrate the accuracy of this procedure. The 


results were as follows: 


NSFO JP-5 AVGAS Tg T DIFFERENCE 
(BBLS) (BBLS) (BBLS) | CALCULATED ACTUAL OBSERVED 
5681 3535 Tee 108 96 re 
4139 12452 0 106 99 +7 


4595 0 0 90 158 -68 
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If this sample can be considered representative, it illustrates quite 
clearly that magnitudes of the errors that might be expected. 

As expected and in agreement with the scatter diagram, the 
standard estimate errors for the cases involving the AE as server were 


less than those involving AOs in four out of six cases. 
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Table IX. 


Results of Linear Regression Analysis on AO and AE 
Server Situations 








AQ-CVA/CVS (Day) 


Sample size 

Percent variance explained 
Multiple correlation coefficient 
Standard error of the estimate 


a - intercept value 
Mean Std. 
Dependent variables (BBLS) Dev. 
NSFO 6577.9 3460.5 
JP-5 6267.8 4638.2 
AVGAS 273.2 098.0 


AO-CVA/CVS (Night) 

Sample size 

Percent variance explained 
Multiple correlation coefficient 


Standard error of the estimate 
a - intercept value 
Mean Std. 
Dependent variables (BBLS) Dev. 
NSFO 6939.5 4579.7 
JP-5 5454.8 3136.9 
AVGAS 258.3 514.4 


AO-CLG/DLG/DDG (Day) 

Sample size 

Percent variance explained 
Standard error of the estimate 


a - intercept value 
Mean Star 
Dependent variable (BBLS) Dev. 
NSFO 2025.1 879.0 


AO-CLG/DLG/DDG (Night) 

Sample size 

Percent variance explained 
Standard error of the estimate 


a - intercept value 
Mean Std. 
Dependent variable (BBLS) Dev. 
NSFO 2111.0 1470.9 


30 


45 
36.4 

604 
27.0 min. 
68.1 min. 


Reg. Std. Err. Partial Corr. 
Coe. Reg. Coe. Coefficients 


0.0050 0.00125 Os6s 
0.0014 0.00093 a 
0.0062 0.00639 aa is 


20 

22.6 
4757 

32.8 min. 

54.9 min. 


Reg. Std. Err. Partial Corr. 
Coe. Reg. Coe. Coefficients 


0.0027 OnO0 Is 0.36 
-0.0013 0.00251 -0.13 
Os02]2 0.01489 0.34 

51 

43.2 

14.4 min. 

24.8 min. 


Reg. Std. Err. Paria l Cove 
Coe. neg. coe. Cocriiciene 


0.0141 0.00230 0.66 
49 
eae 
24.0 min. 
56.5 min. 
Reg. Sta. Evy Partial) Corr. 
Coe. Reg. Coe. Coefficient 
0.0054 0.00236 0.32 


Table IX. (cont.) 


AQ-DD (Day) 
Sample size 99 
Percent variance explained (Ae 
Standard error of the estimate 30.5 mine 
a - intercept value 25. Senin. 
Mean Std. Reg. Std. Err.- PartidileGeny. 
Dependent variable (BBLS) Dev. Coe. Reg. Coe. Coefficient 
NSFO 132099" 533.2" “ORO805 0.00583 0.47 
AQ-DD (Night) 
Sample size 40 
Percent variance explained 0202 
Standard error of the estimate 60.8 min. 
a - intercept value 80.9 min. 
Mean Sid, Reg. Std. Err. Partial €omr. 
Dependent variable (BBLS) Dev. Coe. Reg. Coe. Coefficient 
NSFO 1490.7 449.0 0.0020 0.02164 0.0151 
AE-CVA/CVS (Day) 
Sample size 30 
Percent variance explained 85.5 
Standard error of the estimate 23.4 min. 
a - intercept value -30.4 min. 
Mean , Std. Reg. Std. Err. Partial Corr. 
Dependent variable (S.T.) Dev. Coe. Reg. Coe. Coefficient 
AE Material 207.2 126.9" ‘0.4882 0.03419 OF9ZS 
AE-CVA/CVS (Night 
Sample size 30 
Percent variance explained 6822 
Standard error of the estimate 31.9 min. 
a - intercept value 48.6 min. 
Mean Sid: Reg. Std. Err. Partial Corr. 
Dependent variable (S.T.) Dev. Coe. Reg. Coe. Coefficient 
AE Material 285.6 122.0 0.3%6 0.04861 0.826 


* S.T. denotes Short Tons. 
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Table IX. (cont. ) 


AE-CLG/DLG/DDG (Day) 


Sample size 30 
Percent variance explained 86.7 
Standard error of the estimate 20.7 min. 
a - intercept value 80.3 min. 
Mean Std. Reg. Std.Ghvae. @Partialm@err . 
Dependent variable (S.T.) Dev. Coe. Reg..Coe. Coefficient 
AE Material 50.7 Sone 0.965 0.07144 0.93 


AE-CLG/DLG/DDG (Night) 


Sample size 25 
Percent variance explained 17.8 
Standard error of the estimate 62.5 min. 
a - intercept value 67.4 min. 
Mean Sed: Reg. Std. Err. Partial Corr. 
Dependent variable (S.T.) Dev. Coe. Reg. Coe. Coefficient 
AE Material 90.8 Gone leoo7 ~ %.1489 882 
AE-DD (Day) 
Sample size 30 
Percent variance explained Tons 
Standard error of the estimate 7 oeunin. 
a - intercept value 23.7 Write 
Mean Sed. Reg. Std. Enn. Pantial Gorr. 
Dependent variable (S.T.) Dev. ~. Coe. Reg. Coe. Coefficient 
AE Material 18.0 16.0 0.568 0.0543 0.892 


AE-DD (Night) 


Sample size 30 
Percent variance explained 46.3 
Standard error of the estimate 39.1 min. 
a - intercept value 36.4 min. 
Mean Std. Reg. Std. Err. Partialecorr. 
Dependent variable (S.T.) Dev. Coe. Reg. Coe. Coefficient 
AE Material 2/.9 21.4 limoy/ 0.34 0.68 
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VI. CONCLUSIONS 


It has been shown that the gamma distribution, with parameters 
determined from operational data, can be used to describe the probability 
distribution of CONREP service times for a large number of practical 
situations. With minor adjustments in the gamma parameter estimates, 
the Erlang distribution, with its increased flexibility in analytical 
models employing Laplace transforms, can be used effectively for the 
Same purpose. 

It is conceivable that previous work in analytical and computer 
Simulation studies which assumed exponentially distributed CONREP service 
times can be extended to incorporate this new information, with a result- 
ant increase in the accuracy of these models in depicting actual under- 
way replenishment operations. 

The second part of the analysis indicates that, with few excep- 
tions, estimates of service time utilizing linear regression techniques 
would be too inaccurate, as indicated by the standard error of the 
estimates derived from the data, to serve the purposes for which such 
a model might be used. Although certain cases appeared to be quite 
accurate, the overall results discourage enthusiasm for further effort 
along these lines. It is noted that some cases appear to lend then- 
selves to non-linear estimates, such as might be obtained through poly- 
nomial regression, since some curvature can be detected in several of 
the scatter diagrams. A larger data base than that which was available 
for this study would be necessary to make meaningful conclusions along 


these lines. 


on 
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APPENDIX A 
The Gamma Distribution 


If a process consists of r successive events and if the total 
elapsed time of this process can be regarded as the sum of r independ- 
ent exponential variates each with parameter i, the probability 
distribution of this sum will be a gamma distribution with parameters 
A andr. A form of the gamma distribution may be fitted to many 
positively skewed distributions of statistical data. 

The gamma distribution is described by the following density 
function 

Gl Xcite). = ae 
where A > 0, r > 0, and x are non-negative. 

The sum of r (where r is a positive integer) exponential 
variables each having the same parameter A is also called an Erlang 
distribution. Mathematically, the Erlang distribution is a convolu- 
tion of r exponential variates, i.e., the distribution of the sum of 
r exponential variables. Because of the restriction on r, a process 
which can be represented by an Erlang distribution is easier to work 
with in an analytic sense, especially if Laplace transform techniques 
are utilized. 

The expected value and variance of the gamma distribution are 
given by 

EX 


Wl 


“le ~|3 
NO 


VX = 


From these relationships it follows that 
(EX)? 
VX 


EX 
VX 





Therefore, a gamma distribution, like a normal distribution, 
can be specified by its mean and its variance. 
If r = 1, the gamma distribution is identical to the exponen- 


tial distribution 


As r increases, the gamma distribution approaches a normal 
distribution asymptotically. r is often referred to as the "shape" 
parameter, while A is Known as the "scale" parameter. Figure 2 
illustrates the graphs of g(x;r,A) for several values of r and \ = 1. 
In the gamma distribution the ratio of the standard deviation (the 
square root of the variance) to the mean (expected value) is the 
reciprocal of the square root of r. When this ratio is 1/3 or 
smaller (i.e., when r is 9 or greater), the gamma distribution closely 
approximates a normal distribution with the same mean and standard 
deviation. 

The cumulative distribution function does not exist tn explicit 
form for the gamma distribution, although the values of the so-called 
incomplete gamma function have been tabled by Pearson [7]. The computer 
Subroutine entitled GAMDIS was constructed for and used in this analy- 
Sis. The subroutine employs numerical integration techniques involving 
the trapezoidal rule [8]. GAMDIS has the function of computing proba- 
bilities and expected numbers necessary in making tests of hypotheses 


in chi-square goodness-of-fit tests. Because of the potential usefulness 
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of this subroutine to others engaged in future efforts of this nature, 
a copy of the subroutine, written in FORTRAN with explanatory comment 


cards, is included as Appendix B. 
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APPENDIX B 
Subroutine GAMDIS 


Subroutine GAMDIS is called by specifying the estimates of the 
gamma distribution, R and XLAMDA. These estimates are computed from 
the sample mean and sample variance of the service: times under con- 


Sideration in accordance with the following relationships: 


p= p = (sample mean)* 


sample variance 


) = Sample mean 


XLAMDA Sample variance 


Other calling parameters for the subroutine include NBU, the number of 
units in the sample, and NBCL, the number of cells for which expected 
numbers of units having appropriate service times are to be calculated 
and labeled by vector EXPECT. 

The subroutine performs the following computations: 

1. Defines the limit of the total area for which numerical 
integration will be utilized (LIMIT = NBCL*10). 

2. Defines constants based on given values of R and XLAMDA 
which will be used in specifying the gamma density function (F) used 
in numerical integration. 

3. Sets up subintervals necessary to obtain an accurate 
measurement of the probabilities associated with each interval. 

4. Compute the probabilities mentioned jin 3., (SUM2), plus 
the remainder of probability associated with the last interval (SUM4). 

5. Computes the expected number of units for each cell using 
NBU, SUM2, and SUM4. 

6. Returns the computations of 5. to the main program via 


vector EXPECT. 
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APPENDIX C 


Scatter Diagrams of Operational CONREP Data for Situations 
Involving AO, AE Servers, CVA/CVS, CLG, and DD Customers 
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FIGURE 3. 
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